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In heterogeneous system, the correspondence between calculated and measured quantities, such 
as the conductivity or the resistivity, is not obvious since the former ones are local quantities whereas 
the latter ones are often average values over the sample. In this report, we show explicitly how the 
correspondence can be done in the case of magnetic multilayers. 
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In the linear response regime, the electric current J at 
positioryj- is related to the electric field E at position r' 
througto 
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where i and j refer to the space directions {x, y, z}. The 
conductivity <7ij depends of both r and r' (two-points 
conductivity), it is a local quantity which is expressed 
in the Kubo formalismo as a current-current correlation 
function 

a.,(r,r') = lim - / dt e^* ( [ j(r, t) ; j(r', 0) ] ) , (2) 

where j is the current density and the bracket [ A ; B ] — 
AB — BA. The notation (■ • ■) refers to the configura- 
tional average. The question we address in this report 
is how to link the calculated quantities given by Eq. (||) 
which are local (position dependent) with the measured 
quantities which are mostly non-local. Indeed, the size 
of the electric contacts used to measure the conductivity 
(by measuring the current and the voltage) are gener- 
ally quite large (from /im to mm). Then the measured 
quantity is a kind of average over a part of the sample. 
Strictly speaking, it is rather an effective value than an 
average value. The relation between the measured and 
calculated quantities varies strongly with the geometry 
of the system. 

Let us start with a three-dimensional homogeneous 
system in the stationary regime. In this case, the electric 
current and electric field do not dependent of the posi- 
tion, thus Eq. (j^) reduces to Ji — ^ 
have introduced the effective conductivity 



tjEj where we 



dr' CTy (r,r') 
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which is nothing else than the average value over the sam- 
ple. The r-dependence in the r.s.h. ofEq. (||) is irrelevant 
since the system is invariant by translation. Then, for an 
homogeneous system, the calculated (two-points conduc- 
tivity (Tij{r, r')) and the measured (effective conductivity 
aij) quantities are related by a simple relation. 

It is not more the case for heterogeneous systems be- 
cause the electric current and electric field are not uni- 
form through the sample. In this report, we treat the 



case of a multilayer formed by a superposition of identi- 
cal cells which have the dimension {ao, ag, Luq} where ag 
is the interatomic distance and L the number of layers 
in the cell (see Fig. (la)). We have still an invariance by 
translation in each direction but the periodicity is differ- 
ent in comparison to the three-dimensional homogeneous 
system. The indication r of the position of one site must 
be replaced by the indication R of the position of the 
cell plus the indication I of the position of the site in the 
cell (/ varies from 1 to L). For a multilayer, it is more 
appropriate to write Eq. (|^) under the form 
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(r,r')£;HR'), 
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because we have a translational invariance with respect 
to the vectors R and R', thus the electric current and 
electric field do not depend of the cell position. Follow- 
ing the same procedure than before, we write 



where 



Jih — 
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dR' alf^R, R'). 
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However, we can not design this quantity as the effective 
conductivity since it is a local quantity through the h 
and I2 indices. Since a multilayer is formed by a succes- 
sion of layers made of different material, it is not possible 
to assume that the electric current and electric field are 
layer independent. To get the effective conductivity, we 
have to go further. In the liter|a,t|«re, only approximate 
expressions have been proposeda'B. They are controlled 
by the relative value of two lengths: the mean- free-path 
A and the average thickness d of the uniform layers in 
the cell. When A ^ d (homogeneous limit), the electrons 
can go through many layers without experience any scat- 
tering. They do not feel the difference between the layers 
and thus the cell can be considered as an homogeneous 
system. The layer dependence of the conductivity is not 
relevant and the multilayer can be treated as a three- 
dimensional homogeneous system (see Eq. (^)). When 
X ^ d (local limit), the electrons are scattered many 
times before they go out of one layer: the current in 
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one layer depends only weakly of the electric field in the 
other layers. As a consequence, the contributions of the 
two-points conductivity o^^^^ are negligible, and Eq. (||) 
reduces to J\ = which is still layer dependent. 

Thus, in the local limit, the effective conductivity can not 
be directly obtained. A general derivation of the exact 
effective conductivity valid in any regimes would be more 
appropriate. In this report, we present such a derivation. 
The starting point is Eq. where the two-points con- 
ductivity cr-j'^ links the local electric current J-^ to the 

local electric field E^^ . The conductivity is a (3L x 3L) 
matrix. By inversion of this matrix, we get the matrix 
elements p'j'^ (two-points resistivity) which link the local 

electric field E\'^ to the local electric current j'^ 



(7) 



The important thing now is to consider the experiments, 
in particular the geometry of the system and the way 
how the contacts are made. We distinguish two different 
geometries : the CIP geometry (current in the plane of 
the layers, see Fig. (lb)) and the CPP geometry (current 
perpendicular to the plane of the layers, see Fig. (Ic)). 
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FIG. 1. (a) Schema of the unit cell for a multilayer; 
(b) CIP geometry; (c) CPP geometry. 

In CPP geometry, the current is perpendicular to the 
layers. Due to current conservation, the perpendicular 
current is uniform through the cell (i.e., j' — Jz) and 
the measured value is Jz. In CIP geometry, the current 
is parallel to the layers. Since the layers are different, 
the in- plane currents and Jy are layer dependent. As 
it is shown in Fig (Ic), the contacts cover several lay- 
ers. We can assume that they cover at least a whole 
cell and that they do not have any influence on it, then 
what is measured are the average currents Jx = X); •^x/-^ 
and Jy = 'YliiJy/L- For both geometries, the electrical 
field can have parallel and perpendicular components. 
Since the layers are different, the three components of 
the electric field should be layer dependent. However, the 
Maxwell equation in the stationary regime (V x E = 0) 
combined with the fact that the electric field compo- 
nent E^^ is uniform in the xy-plane impose that in-plane 
components of E are layer independent (i?^ — E^ and 
Ey = Ey). Thus, Ex and Ey correspond to the mea- 



sured quantities. On the contrary, E\ stays layer depen- 
dent and what is measured is the average electric field 

By means of some transformations on Eqs. (|5|) and ([^), 
we express J^, Jy and E\ with the help of the uniform 
quantities Ex^ Ey and Jz- For the moment, all the indices 
are kept. Eq. (||) can be written as 



^From Eq. (0), we have 
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These three equations can be written under the matrix 
form 




(11) 



where / is the (L x L) identity matrix and where we have 
introduced the following notations: the L-components 
vector Ai {A = J or E) and the (L x L) matrix 
(a = cr or p) define as 



A,, 



From Eq 
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We inverse the first matrix in the r.s.h (it can be done 
numerically) and perform the multiplication of the two 
matrices, thus we get 
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where Sij and are (-L x _L) matrices defined through By using the fact that E^, Ey and are layer- 
independent, we can sum over the columns and reduce 
^xy \ / pxx Pxy ^ \ " the initial size (3L x 3L) of the matrix which appears in 

^yx Syy = Pyx Pyj^ ] . (15) Eq. (|i|) to thc sizc (3i X 3) 
p,, / \ 







'j'x 




n 

Jy 




Uy 




El 









Ah 



E 
E 



^xx 
Ah 



E 

E/2,;3 P 



„Lh 

h 



\-Y.h,hPzz'^\ 



zx 



El 



Eia ^xy 
Ah 

■ih yy 



El. 



E 



E 



^Lh 

h ^yy 



E 



h,i. 



zy 



rfLh I3I2 
I zz zy 



Eh 



Ah J^ih 



l'3'2 _1_ AI3 li'i^'i \ \ 
"XX rxz ' xy ryz } 



Eijjis ^^xxPxz' + ^xyP. 

E {Ah Jih I Ah J3I 
h,h \ yx Pxz ^ ^yy Hyz 



Lh /^hh 
xy rxy 



E 



h,h \ y^ 



Ih 



EhP- 



Eh Pzz^ 



^Lh /^hh 

yy i^yz 




J 



(16) 



As we are interested by the average currents = Ei ^x/-^j Jy = Ei Jy/-^ ^^^'^ the average electric field = EiEi/L 
(it is what we can get experimentally in such a system) , we reduce the (3L x 3) matrix which appears in Eq. dlq) to 
a (3 X 3) matrix 




I ^xx 1 1 

WPzz^zx I 
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where we have introduced the definition 
equations can be written under the form 



Eu, 



/L in order to simplify the notations. This system of 
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From Eq. (18), the effective conductivity tensor a is 
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The inversion of the first matrix and the product of the 
two matrices lead to the following expressions of the ma- 
trix elements for the effective conductivity 
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The inversion of the first matrix and the product of the 
two matrices lead to the fohowing expressions of the ma- 
trix elements for the effective resistivity 
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where we have introduced the denominator D — 
I|s2;a;||||syy|| - ||s^y||||sya:||- Eqs. (|2^) and (|2|) are valid 
in the general case. We shall now give reduced expres- 
sions in some particular cases. 

Without magnetization and for a cubic lattice, the 
symmetry imposes cr-^'j = p'^^ =0. In addition, we 

have s^i'^ = al]'^^ for i e {x,y} and p^ij^ = p^ij^ (see 
Eq. dlq)). Thus, the effective conductivity and resistiv- 
ity tensors (Eqs. ( p^ and (pl|), respectively) reduce to 
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When the magnetization is along the x-direction, we 
have a^^ = a"^ = p'^^^ = = for * e {y,z} 

(and similarly for s and p). They are no particular rela- 
tions between s'j'^ and trl]'^ or between p'j'^ and p'^^ 
Thus, the effective conductivity and resistivity tensors 
are equals to 
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where D = ||sa;a;|| ||si/j/||. The results with a magnetiza- 
tion along the y-direction can be obtained from Eqs. ( |25| ) 
and (|2^) by exchanging x and y indices. 

When the magnetization is along the z-direction, we 
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= P.. = Pzt = for z e {x,y}. In 
for {ij} e {x,y} and 
(see Eq. (|l5|)). Thus, the effective conduc- 



have 

addition, we have sl)''^ — al]^ 



tivity and resistivity tensors are equals to 




(27) 
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where D = HcTxi^II Ikyyll " ||o-2;y|| Iky2;||- The fact that for 
CIP, the effective conductivity a^j (where {i,j} G {x,y}) 

is simply equal to ||cry || — I2 '^Ij''^^ / ^ ^^^^ 
for CPP, the effective conductivity a^z is equal to 
VIIp^^II = Pi.'V-^) has been widely used in 

the literatureoQiJ. Assuming the local limit, it is thus 
possible to model the multilayer as a network of resistors 
in series (in the-JDIP geometry) or in parallels (in the 
CPP geometry )0l3. 
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